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Ionospheric  Scintillation  Effects  on 
VHF-UHF  Communication  Systems 


I.  INTRODUCTION 

Satellite  communication  links  at  UHF  can  be  subject  to  the  effects  of  ionospheric 
scintillations,  which  are  principally  related  to  the  occurrence  of  F-layer  irregulari- 
ties. Scintillations  cause  both  enhancements  and  fading  about  the  median  level  as 
the  radio  signal  transits  the  disturbed  ionospheric  region.  When  scintillations  occur 
which  exceed  the  fade  margin,  performance  of  the  communications  link  will  be  de- 
graded. This  degradation  is  most  serious  for  propagation  paths  which  transit  the 
auroral  and  equatorial  ionospheres.  The  degree  of  degradation  will  depend  on  how 
far  the  signal  fades  below  the  margin,  the  duration  of  the  fade,  the  type  of  modula- 
tion and  the  criteria  for  acceptability. 

The  occurrence  of  scintillations  has  been  studied  for  several  years  and  its 

morphology  has  been  documented  for  the  auroral,  mid-latitude  and  equatorial 
12  3 

regions.  ' * This  report  concentrates  on  the  analysis  of  specific  periods  of 
intense  scintillations  and  applies  the  results  to  the  evaluation  of  communication 
system  performance. 

(Received  for  publication  3 May  1976) 

1.  Aarons,  J.,  Whitney,  H.  E, , and  Allen,  R.  S.  ( 197  1)  Global  morphology  of 

ionospheric  scintillations,  Proc.  IEEE  ^9:159, 

2.  Crane,  R.K,  (1974)  Morphology  of  Ionospheric  Scintillation,  Technical  Note. 

Lincoln  Laboratory,  1974-29,  Air  Force  Contract  F19628-73-C-0002. 

3.  Aarons,  J.  (1975)  Global  Morphology  of  Ionospheric  Scintillations  II,  AFCRL- 

TR-75-0135. 
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2.  DATA  - FORMAT  AND  ANALYSIS 


|.)nta  from  Huancnyo,  Peru  ami  Nurssarssunq,  Greenland  were  used  for  this 

4 

study.  The  data  were  analyzed  in  15-min  intervals  as  earlier  studies  had  shown 
that  this  length  was  a good  choice  for  scintillation  studies  using  synchronous  satel 
lites.  The  periods  covered  by  each  set  of  data  are  listed  in  Table  1. 


Table  !.  Summary  of  Data  Hnse 


Site 

Freq  (MHz)  • 

Date 

Time  (UT) 

No.  of 
15-min. 
Intervals 

Peru 

137/360 

2 Nov  1974 

0315  to  0500 

7 

Peru 

137/360 

5 Nov  1974 

0110  to  0433 

10 

Narssarssunq 

137 

90  .lul  1974 

0500  to  0830 

14 

This  selection  of  data  provided  a total  of  31  periods  of  quite  intense  scintillations. 
These  data  were  not  selected  to  represent  'worst'  case  conditions  but  only  u sampl- 
ing of  strong  scintillations  that  can  occur  in  the  polar,  auroral  and  equatorial  re- 
gions. The  360  MHz  data  from  Peru  was  at  a frequency  of  direct  interest  to 
AFSATCOM,  while  only  137  MHz  data  was  available  on  magnetic  tape  from 
Nar.ssarssuaq.  The  Peru  data  was  emphasized  because  it  exhibited  the  strongest 
scintillation  and  had  two  frequencies  which  provided  a measure  of  the  spectral  index 
or  frequency  dependence  of  scintillation  for  this  frequency  range. 

The  signal  strength  level  from  satellite  beacons  was  recorded  on  FM  analog 
tape  and  then  digitized  at  six  samples  per  sec  for  computer  processing.  A relative 
calibration  of  the  complete  receiving  and  recording  equipment  was  accomplished 
near  the  time  of  each  scintillation  record. 

A computer  program  was  written  that  provided  the  following  outputs: 

(a)  HIM  plot  of  the  signal  level, 

(b)  cumulative  probability  amplitude  distribution,  edf, 

(c)  fade  rate  distributions, 

(d)  message  reliability, 

(e)  power  spectrum, 

(f)  autocorrelation,  and 

(g)  cross-correlation. 

4.  Whitney,  H.E.,  Aarons,  J. , Allen,  R.S.,  and  Seemann,  D.  R.  (1972)  Estima- 
tion of  the  cumulative  amplitude  probability  distribution  function  of  ionos- 
pheric scintillations.  Radio  Science,  7,  No.  12:1095-1104. 
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Aii  example  of  each  of  the  outputs  Is  shown  in  Figures  1 through  8 for  both  the 
137  Mil*  and  360  Milt  data, 

A IS-min  section  of  the  chart  recording  at  each  frequency  and  the  resulting  cdf 
is  shown  in  Figure  1,  The  cdf  Is  a first  order  statistic  and  is  useful  for  defining 
the  minimum  margin  requirements  for  the  communications  link  of  nondiversity 

5 

systems.  The  various  scintillation  indices  defined  by  Briggs  and  Parkins,  that  Is, 
S j,  S,„  S and  S4  are  calculated  for  each  15-mln  samp]*  s which  depends  on  the 
rms  value  of  the  power,  is  the  most  commonly  used  Index,  The  Nakagamt  m- 
distribution  has  been  shown  to  be  practically  useful  for  describing  the  effects  of 

7 

scintillations  on  satellite  communication  links.  Christopher  has  described  a 

useful  application  of  the  moments  of  a log-Nakagaml  distribution.  Also  calculated 

2 

for  each  15-mln  interval  are  the  Nakngami  m-value,  (m  --=  1/S^);  a Chi-square  fit- 
ness test  of  the  experimental  data  to  both  Nakagaml  and  log-normal  distributions 
over  the  2 to  96  percentile  range;  and  the  first  four  moments.  MU,  of  the  log- 
Nakagami  pdf.  The  circles  in  Figure  1 are  points  for  a Rayleigh  distribution  (m=  1) 
and  show  the  agreement  with  the  experimental  datR. 

In  addition  to  the  information  on  the  r nplitude  of  the  fades  which  is  given  by  the 
cdf,  statistical  descriptions  of  the  fading  rate  is  needed  in  order  to  fully  character- 
ize the  effects  of  scintillations  on  the  communications  channel,  Information  on  the 
fading  rate  can  either  be  produced  by  level-crossing  techniques  or  by  Fourier  tech- 
niques which  produce  the  power  spectra  and  time  correlation  functions.  The  level- 
crossing data  give  measurements  of  the  distribution  of  the  fades  and  enhancements 
at  various  levels  relative  to  the  median  level.  The  duration  and  separation  of  the 
fades  are  counted  and  tabulated  against  time  intervals  for  1 dB  changes  in  signal 
level  down  to  -10  dB.  Intervals  shortei  than  0.  25  sec  sec  and  longer  than  60  sec 
were  not  resolved.  Figure  2 shows  the  distribution  of  fades  below  several  levels 
for  both  137  MHz  and  360  MHz,  For  this  particular  sample  the  occurrence  of  fades 
of  a given  duration  at  137  MHz  was  approximately  twice  that  occurring  at  360  MHz, 
The  presentation  of  enhancements  of  signals  for  various  levels  can  be  done  in  a 
similar  manner. 

Figure  3 shows  a plot  of  the  message  -eliability  for  both  frequencies.  This 
given  an  estimate  ot  the  increase  in  margin  which  is  required  over  the  value  speci- 
fied by  the  cdf  to  obtain  a given  probability  of  receiving  perfect  messages.  Message 
reliability  waa  measured  by  determining  the  number  of  time  intervals  or  messages 

5.  Briggs,  B.H. , and  Parkins,  I,  A.  (1963)  On  the  variation  of  radio  star  and 

satellite  scintillations  with  zenith  angle,  J,  Atmos.  Terr,  Phy.  , 25:330.  366. 

6.  Nakagami,  M.  (1960)  Statistical  methods  in  Radio  Wave  Propagation,  edited 

by  W,  C.  Hoffman,  pp.  3-36,  Pergamon,  New  York 

7.  Christopher,  P.  (197  5)  Moments  of  a Log-Nakagami  Distribution  and  Applica- 

tion, Mitre  Corporation  Memo  D9I-M-40<24, 
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the  fading  rate  increases  the  curves  will  fall  off  more  rapidly.  If  the  fading  rate 
becomes  so  rapid  that  the  inter-fade  separations  are  less  than  the  message  length, 
no  messages  will  be  received  without  error. 

Figure  4 shows  the  137  MHz  and  360  MHz  power  spectra  for  this  data  sample. 
The  rate  of  the  scintillation  fading  is  characterized  by  the  power  spectrum.  In 
general,  the  spectral  shape  consists  of  a relatively  flat  low  frequency  spectrum 
and  a high  frequency  roll-off  with  a slope  of  approximately  f 1 for  intense  scintilla- 
tions. A 'cut-off  frequency  denotes  the  break  point  between  the  low  and  high  fre- 
quency slopes.  Some  spectra  had  a fairly  abrupt  transition  in  the  0.  1 to  1 Hz  range 
between  the  low  and  high  frequency  bands  and  the  Intersection  of  straight  line  fits  to 
the  slopes  give  the  'cut-off'  frequency.  The  results  are  tabulated  in  Table  2.  Other 
spectra  did  not  exhibit  a 'cut-off'  frequency  and  only  a single  slope  was  measured. 
The  strong  discrete  lines  are  associated  with  the  satellite  spin  modulation. 

Figure  5 shows  the  137  MHz  and  360  MHz  autocorrelation  functions  for  the  data 
sample.  The  autocorrelation  function  is  another  way  of  characterizing  the  rate  of 
scintillation  fading.  It  is  the  Fourier  transform  of  the  power  spectrum  and  there- 
fore has  a width  which  is  inversely  proportional  to  the  bandwidth  of  the  power 
spectrum  where  the  bandwidth  is  defined  by  the  'cut-off'  frequency.  The  correla- 
tion interval,  r is  necessary  to  evaluate  the  effectiveness  of  time  diversity.  Coding 
techniques  can  be  an  effective  means  of  achieving  time  diversity  improvement.  The 
correlation  for  the  full  period  is  plotted  and  in  addition  the  first  16  sec  is  expanded. 
The  correlation  interval,  r , is  printed  out  for  coefficients  of  approximately  0.  5, 

0 and  at  the  first  inflection  point. 

Figure  6 shows  the  cross-correlation  function  of  the  137  MHz  and  360  MHz 
channels.  It  indicates  the  effectiveness  of  applying  frequency  diversity  techniques 
to  overcome  the  degradation  due  to  scintillations.  The  cross-correlation  coefficient 
was  less  than  0.  2 for  this  sample  and  also  was  low  (<.  4)  for  two  other  data  samples 
that  were  analyzed.  It  is  expected  that  the  cross- correlation  values  would  be  low 
for  this  wide  a frequency  difference  (137/360  MHz)  under  conditions  of  intense 
scintillations.  Since  Schwartz  et  al  indicates  only  a very  small  improvement  in 
diversity  gain  for  correlation  coefficients  between  0.  5 and  0 the  computation  of  the 
cross-correlation  function  was  not  done  for  all  the  two  frequency  data.  From  the 
standpoint  of  the  needs  of  the  AFSATCOM  program  it  would  have  been  more  desir- 
able ,o  have  been  able  to  measure  the  cross-correlation  interval  for  a smaller  fre- 
quency separation.  However,  the  137/360  MHz  data  were  the  only  available  two 
frequency  data. 


8.  Schwartz,  M. , Bennet,  W.  R. , and  Stein,  S.  (1966)  Communication  Systems  and 
Techniques,  McGraw-Hill.  ~ " 
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Figure  6.  Typical  Crosscorrelation  Data  for  Intense  Scintillations  at  137  MHz  an< 


3.  RESULTS 


Table  2 is  a summary  of  the  Important  statistics  of  the  analysis  of  the  31  data 
samples.  The  table  lists  the  Nakagami  m-parameter;  the  correlation  interval,  T, 
for  an  autocorrelation  coefficient  p = 0.  5;  the  'cut-off'  frequency  and  the  slopes  of 
the  power  spectrum;  the  calculated  value  of  the  spectral  index  for  samples  having 
two  frequencies;  and  the  calculated  value  of  the  product  of  the  'cut-off  frequency 
and  the  correlation  interval,  r. 

The  Nakagami  m-parameter  is  related  to  the  rms  value  of  the  intensity  of  the 
scintillations.  It  completely  describes  the  distribution.  Reference  4 has  shown 
that  the  distribution  function  of  ionospheric  scintillations  closely  approximates  the 
theoretical  Nakagami  distribution.  While  m can  have  any  value  > 0.  5,  all  the  data 
that  has  been  analyzed  Indicates  that  m = 1 is  the  limiting  case  which  corresponds  to 
the  Rayleigh  distribution.  The  circles  in  Figure  1 are  points  for  a Rayleigh  distrib- 
ution, (m=  1),  and  show  the  agreement  of  the  137  and  360  MHz  data  with  the  Naka- 
gami distribution.  The  Chi-square  test  indicates  that  there  is  a better  fit  of  the 
data  to  the  Nakagami  distribution  than  to  the  log-normal  distribution.  Figure  7 is 
a plot  of  the  Nakagami  distribution  in  cumulative  form  for  several  values  of  m.  The 
experimental  data  indicates  that  the  Rayleigh  distribution  or  the  m=  1 curve  in  Fig- 
ure 7 is  a limiting  case  for  scintillations  in  the  200  to  400  MHz  band.  Figure  7 can 
therefore  be  used  to  estimate  the  expected  fading  under  'worst'  case  conditions  for 
various  percentiles. 

Reference  4 has  shown  that  the  spectral  index  or  frequency  dependence  of  scin- 
tillations can  be  related  to  the  m-parameter  in  the  following  way; 

m, 

log 

V 


m 


m, 

4 

T T 


log 

2 


where  n is  the  spectral  index  and  m.  and  m„  are  the  Nakagami  values  measured 
at  the  two  frequencies,  f.  This  assumes  that  the  scintillation  index  or  m remains  a 
constant  power  of  frequency  over  the  measurement  range.  Table  2 lists  the  calcu- 
lated values  of  the  spectral  index  for  the  17  sets  of  two  frequency  data  and  the  dis- 
tribution is  plotted  in  Figure  8.  While  the  scintillations  at  360  MHz  did  not  have  the 
intensity  to  actually  measure  a value  of  r Jm  = 0,  it  is  expected  that  this  could  be 

achieved  as  the  'worst'  case  condition.  There  were  several  measured  values  of  v 

'm 

in  the  lowest  range  of  0.  2 which  indicates  very  little  decrease  of  scintillation 
intensity  with  frequency.  This  is  another  indication  that  the  'worst'  case  fading  in 
the  200  to  400  MHz  band  can  be  approximated  by  a Rayleigh  distribution. 


PERCENT  OE  TIME  SIGNAL  INTENSITY,  X,  EXCEEDED  ORDiNATE 


Figure  7.  Cumulative  Form  of  Nakagami 
m- Distribution:  Valid  for  m>  0.  5 


In  addition  to  the  information  on  the  depth  of  fading,  information  on  the  rate  of 
fading  i3  also  required  to  fully  understand  the  effects  of  scintillations  on  the  com- 
munications channel.  The  statistics  on  the  fading  rates  has  been  measured  and 
presented  in  three  ways;  (1)  the  autocorrelation  function,  (2)  the  power  spectrum, 
and  (3)  the  tabulation  of  the  results  of  the  level  crossing  technique  in  the  form  of 
the  duration  of  fades  below  and  increases  above  given  levels  and  message  reliability. 
Since  time  diversity  is  effectively  achieved  in  the  application  of  coding  techniques 
to  the  AFSATCOM  system,  it  is  felt  that  the  autocorrelation  information  is  the 
easiest  to  apply  In  evaluating  the  fade  rate  data. 
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SPECTRAL  INDEX, 


Figure  8.  Occurrence  of  Spectral  Index  (tj  ) 
Values  From  a Data  Base  of  17  Pair  of  15-min 
Amplitude  Distributions  of  137  and  360  MHz 
Scintillations;  m = Nakagami's  Distribution 
Parameter 


When  correlation  data  is  available  it  can  be  used  to  determine  the  improvement 
in  performance  that  can  be  obtained  through  the  application  of  diversity  techniques. 
Autor' rrelation  data  can  be  used  to  evaluate  time  diversity  techniques;  cross- 
correlation  data  from  multifrequency  measurements  can  be  used  to  evaluate  the 
effectiveness  of  frequency  diversity;  and,  if  spaced  receiver  measurements  are 
available,  then  the  cross-correlation  information  can  evaluate  space  diversity. 
Figure  9 is  adapted  from  Reference  8 and  shows  the  improvement  in  performance 
that  is  possible  with  dual-diversity  techniques.  The  probability  of  achieving  diver- 
sity gain  is  given  for  several  values  of  p,  the  correlation  coefficient.  It  is  based 
on  slow,  multiplicative  Rayleigh  fading  and  equal  signal-to-noise  ratios  in  both 
branches  of  a dual  diversity  system.  Reference  8 used  I pi , the  magnitude  of  the 
complex  cross-correlation,  while  Figure  9 uses  the  envelope  cross- correlation. 
Most  of  the  diversity  improvement  is  achieved  by  the  time  p = 0.  6.  For  example 
the  improvement  at  the  1 percent  point  is  approximately  8 dB  for  p = 0.  6 and  10  dB 
for  p = 0 or  complete  decorrelation. 

The  correlation  intervals  for  p = 0.  5 of  the  17  sets  of  137/360  MHz  data  are 
plotted  in  Figure  10  to  show  the  relationship  that  holds  for  this  very  limited  set  of 
data.  For  small  delays  the  360  MHz  value  is  approximately  twice  the  137  value. 

At  longer  delays  the  values  tend  to  be  the  same. 
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The  variation  of  the  correlation  interval  with  Intensity  of  scintillation  Is  shown 
In  Figure  11,  The  measured  values  at  each  frequency  for  the  17  sets  of  data  are 
plotted  against  the  Nakagaml  m-parametcr  which  was  measured  from  the  corres- 
ponding cdf.  An  ellipse  has  been  drawn  to  encompass  the  data  at  each  frequency 
to  give  the  apparent  variation  with  m. 


Figure  11.  Variation  of  the 
Correlation  Interval  With  Intensity 
of  Scintillations  for  Two  Frequency 
Data  Samples 


The  values  of  the  correlation  interval  for  the  Narssarssuaq  data  were  generally 
greater  than  for  the  Peru  data,  even  though  they  were  both  measured  at  the  same 
frequency,  137  MHz.  It  is  felt  that  this  is  not  a difference  between  the  equatorial 
and  auroral  regions  but  rather  it  is  caused  by  the  Narssarssuaq  data  having  less 
intensity  of  scintillations  than  the  Peru  data. 


For  those  periods  given  in  Table  2 for  which  a cutoff  frequency  could  be  mea- 
sured in  the  power  spectrum,  a constant  was  calculated  from  the  product  of  the 
cutoff  frequency  and  the  autocorrelation  interval,  r , The  values  are  tabulated  in 
Table  2 under  the  heading  and  Kg6Q.  The  averages  of  the  measurements  for 
both  frequencies  were  almost  equal,  Kjg7  = ,37  and  Kggg  = , 35. 

4 CONCLUSIONS  AND  RECOMMENDATIONS 

From  a large  collection  of  observations  a limited  amount  of  VHF/UHF  data  has 
been  analyzed  to  determine  the  amplitude  and  rate  characteristics  of  Intense  scintil- 
lations. The  results  are  in  the  form  of  amplitude  distributions,  level  crossing 
tabulations,  power  spectra  and  correlation  data.  While  a Rayleigh  distribution  was 
not  exactly  measured,  it  was  judged  that  this  would  be  the  limiting  case  up  to  at 
least  400  MIIz.  This  is  in  agreement  with  the  spectral  index  data  which  indicated 
that  fading  of  the  most  intense  samples  was  almost  independent  of  frequency.  The 
autocorrelation  and  power  spectra  data  defined  the  fading  rates  and  provided  a base 
for  the  evaluation  of  time  diversity  techni  es.  Only  three  samples  of  cross- 
correlation (137/360  MHz)  were  process  and  coefficients  less  than  0.4  were 
obtained.  However,  this  frequency  separation  was  really  too  great  to  determine 
the  applicability  of  frequency  diversity  techniques  within  the  bandwidth  normally 
available  to  satellite  communication  systems. 

Much  remains  to  be  done  to  characterize  the  Intense  scintillations  that  are 
obseiwed  in  the  equatorial  region.  Multifrequency  observations  are  necessary  to 
determine  the  frequency  dependence  and  the  evaluation  of  diversity  techniques  in 
the  UHF  range. 
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